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Abstract 



We analyse the relic neutralino density in supersymmetric models with an interme- 
diate unification scale. In particular, we present concrete cosmological scenarios where 
the reheating temperature is as small as 0{1 — 1000 MeV). When this temperature is 
associated to the decay of moduli fields producing neutralinos, we show that the relic 
abundance increases considerably with respect to the standard thermal production. Thus 
the neutralino becomes a good dark matter candidate with 0.1 < Qh'^ < 0.3, even for re- 
gions of the parameter space where large neutralino-nucleon cross sections, compatible 
with current dark matter experiments, are present. This is obtained for intermediate 
scales Mj ~ 10^^ — 10^^ GeV, and moduli masses ~ 100 — 1000 GeV. On the other 
hand, when the above temperature is associated to the decay of an infiaton field, the relic 
abundance is too small. 



1 Introduction 



As it is well known, the lightest neutralino, Xi? is a weakly interacting massive particle 
(WIMP), and therefore a very interesting candidate for dark matter in the universe. In 
fact, many experimental efforts are being carried out in order to detect WIMPs through 
elastic scattering with nuclei in a detector . In this sense the theoretical analysis of the 
neutralino-nucleus cross section o"^o_^ is very important. In particular, these analyses 
in the context of the minimal supersymmetric standard model (MSSM) are usually per- 
formed assuming the unification scale Mqut ~ 10^^ GeV for the running of the universal 
soft supersymmetry (SUSY)-breaking terms. However, it was pointed out recently that 
this cross section is very sensitive to the variation of the unification scale. For instance, 
by taking an intermediate unification scale Mj k, 10^°~^^ GeV the cross section increases 
substantially, being compatible for large regions of the parameter space of the MSSM 
with the sensitivity of current dark matter experiments o"^o_^ k, 10~^-10~^ GeV~^, for 
tan /3 > 3 and m^o ^ 100 GeV. For larger values of the scale, as e.g. M/ = 10^^ GeV, 
a similar result is obtained for tan/5 > 10. Explicit scenarios with intermediate scales, 
arising in D-brane constructions from type I strings, were analysed in ref. 0. Although 
compatibility with the experiments may also be obtained within the usual MSSM scenario 
with the scale ~ 10^^ GeV, it requires large values of tan [3 (tan [3 > 20) or a specific 

non-universal structure of the soft terms 0, ^ ^ §. 

In all the above works the relic neutralino density was also discussed. In these scenarios 
with a large cross section in some regions of the parameter space, generically ^x^h'^ ^ 0.01. 
Of course, this might be a potential problem for the consistency of those regions given 
the observational boundsQO.l < ^xoh"^ ^ 0.3. 

This result is obtained because in the usual early-universe model thermal production 
of neutralinos gives rise to Q^oh'^ oc l/(cr-o'^f), where a-o"^ is the cross section for anni- 
hilation of a pair of neutralinos, v is the relative velocity between the two neutralinos, 
and (..) denotes thermal averaging. Therefore, in this scheme the relic density is inversely 
proportional to the annihilation cross section. Let us recall that crossing arguments, when 
the main annihilation channel is into quarks, ensure that the cross sections of annihilation 
and scattering with nucleons are similar. Thus a large scattering cross section cr^o_j^ leads 
generically to a large annihilation cross section c^o"^, and as a consequence to a small relic 
density. 

However, it is important to remark that this result depends on assumptions about the 
evolution of the early universe. In principle, different cosmological scenarios might give 
rise to different results. To address this question is precisely the aim of this paper. We 
will study the relic density in the context of some non-standard cosmological scenarios. 
In particular, we will show that, when intermediate scales are present, results different 
from the usual ones summarised above may be produced. This is because a low reheating 
temperature, below the freeze-out temperature, can be obtained. We will see that, in the 
case of one of the scenarios, values of the relic density within the observational bounds 

^It is worth noticing, however, that more conservative lower bounds, fl^oh^ w 0.01, have also been 
quoted in the literature. For a brief discussion on this issue see e.g. ref. and references therein. 
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are possible, even for regions of the parameter space with a large neutralino-nucleus cross 
section a^o-iv ~ 10"'' GeV"^. 

The content of the paper is as follows. In Section 2 we will briefly review the usual 
cosmological scenario where thermal production of neutralinos is assumed. Several well- 
known formulas will be explicitly written since we will use them in the discussions of 
the next sections. Then, in Section 3, we will discuss the modifications introduced in 
the relic density analysis by considering non-standard cosmological scenarios in the case 
of intermediate scales. In particular, we will study the situation when an inflation or 
a modulus field produce low reheating temperatures, close to the nucleosynthesis one. 
Finally, the conclusions are left for Section 4. 



2 Thermal production of neutralinos 

Let us briefly review the standard computation of the cosmological abundance of neutrali- 
nos 1^. Neutralinos were in thermal equilibrium with the standard model particles in the 
early universe, and decoupled when they were non-relativistic. The process was the follow- 
ing. When the temperature T of the universe was larger than the mass of the neutralino, 
the number density of neutralinos and photons was roughly the same, n^o oc T^, and the 
neutralino was annihilating with its own antiparticle into lighter particles and vice versa. 
However, shortly after the temperature dropped below the mass of the neutralino, fn^o, 

its number density dropped exponentially, n% cx e ^A^'^ ^ because only a small fraction 

Xi 

of the light particles mentioned above had sufficient kinetic energy to create neutralinos. 
As a consequence, the neutralino annihilation rate F^o = {a'io"'v)n^o dropped below the 
expansion rate of the universe, F^o < H, where H is the Hubble expansion rate. At this 
point neutralinos came away, they could not annihilate, and their density is the same 
since then. This can be obtained using the Boltzmann equation: 

^ + 3Hn,o = -(a|r.) [(n,o)2 - (n-)^] . (1) 

One can discuss qualitatively the solution using the freeze-out condition F^o = (cr|o"f )^?t,^o 

H. Then f^^^o/i^ = {p^o/pc)h'^, where p^o is the current neutralino mass density and Pc is 
the critical density, turns out to be 

^ ^x; H _ m,o (45/67rv/l0) 

(27rV45) g,{Tp) T| {a^rv)F pjso Mp gl'\Tp) Tp {a'^^v)p Pc/so ' ^ ' 



where g^{T) is the effective number of degrees of freedom at temperature T (e.g. including 
all the standard-model degrees of freedom one gets = 106.75), Tp is the freeze-out 
temperature, and sq is the current entropy density. In the second expression we have 
used the fact that H = {7T^g^{T)/90f^^ T^Mp\ with Mp = Mpianck/VS^ ^ 2.4 x 10^^ 
GeV the reduced Planck mass. Taking into account the current value Pc/so — 3.6 x 10~^ 
GeVx/i^, and the typical freeze-out temperature Tp ~ m^o/20, one can write the above 
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expression as 



1 GeV-'^\ ( 100 ^^'^ 



Q^,h' ~ 1.7 X 10"^" — — . (3) 



Since neutralinos freeze-out at Tp ~ m^o/20 << "ni^o, they are non-relativistic and 
therefore the averaged annihilation cross section can be expanded as follows: 

(a^rt;) = as + ap{v') , (4) 

where as describes the s-wave annihilation and ap describes both s- and p- wave annihi- 
lation. Then, eq. (0) can alternatively be written as: 

Q^oh^ ~ 8.8 X 10-"— ^^^^ ^ (5) 

' gj (Tp) [as/xp + Mp/xp^ 

where xp = m^/Tp. 

As it is well known, in most of the parameter space of the MSSM the neutralino is 
mainly pure bino, and as a consequence it will mainly annihilate into lepton pairs through 
t-channel exchange of right-handed sleptons. The p-wave dominant cross section is given 



by m 11 



{a^-v) ^ 8W^-^^— ^ {v^) , (6) 



<1 



where xj^ = m- /tu'^o and a' is the coupling constant for the U{1)y interaction. Taking 
~ m^o ~ 100 GeV, (cr|o"f ) in eq. (P) becomes of the order of 10~^ GeV"^ or smaller. 
Using eq. (|) an interesting relic abundance, ^^oh? > 0.1, is obtained. 

However, in the special regions mentioned in the Introduction, with non-universality 
and/or large tan/3, the lightest neutralino may have an important Higgsino component, 
producing a larger cross section. This is also the case of scenarios with an intermediate 
unification scale. The upper bound for the annihilation cross section, obtained when the 
neutralino is Higgsino-like, is given by [|T^] 

^ ^ " 2 m|o (2 - x^f ' 

where xw = t^wI^\o and a2 is the coupling constant for the SU{2)l interaction. Here 
one is considering that the Higgsino dominantly annihilates into W-boson pairs. Since 
now {a'^^v) in eq. (|^) is of the order of 10"*^ GeV~^, the relic abundance given by eq. 

turns out to be small, ^xoh"^ ~ 0.01, as expected. 
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3 Non-standard cosmological scenarios 



In the standard computation reviewed in the previous section, one is tacitly assuming 
that the radiation-dominated era is the result of a reheat process in the early universe, 
where the reheating temperature Trh is very large, in particular Trh » Tp ~ 10 GeV. 
The scalar field 0, whose decay leads to reheating, is usually assumed to be the infiaton 
field. One can estimate the reheating temperature as a function of the decay width as 

m 

Trh={ , X \v,MpflK (8) 

However, the only constraint on the reheating temperature is Tpu > 1 MeV in order 
not to affect the successful predictions of big-bang nucleosynthesis. This allows in princi- 
ple to consider cosmological scenarios with a low reheating temperature []1T|, T^h < Tp- 
On the other hand, the reheating process can also be associated with the decay of moduli 
fields, as e.g. those appearing in string theory. Thus the relic abundance could receive 



contributions from this source 13 



In what follows, we will show that scenarios with intermediate unification scales are 
explicit examples for the two non-standard cosmological possibilities mentioned above, 
namely, decay of the infiaton and moduli fields producing a low reheating temperature. 
For this analysis eq. (Q) is still valid using F,^ as given by the corresponding scenario. This 
is also true for relation VL^oh? oc l/(cr|o°f) in the case of neutralino production through 
modulus decay. Notice that in this case i^^o oc 1/Tp, as shown in eq. (0), and therefore if 
this mechanism produces a temperature smaller than the typical Tp ^ m^o/20, the value 
of the relic neutralino density will be increased |. We will see below that temperatures as 
required can be obtained in scenarios with intermediate scales. In ref. [|l^ this mechanism 
was applied in order to obtain reasonable values of the relic wino density in anomaly- 
mediated SUSY breaking scenarios, using ~ 100 TeV. In our case, standard masses 
in supergravity scenarios, ^ 1 TeV, will be used. 

On the other hand, in the scenario where the low reheating temperature is obtained 
through an infiaton field, the result for the relic density is quite different from the usual 
one with large Trh- In fact, in certain cases, the usual relation VL^oh? oc 1/ (a^^v) is not 
even valid, and the relic abundance may well be proportional to the annihilation cross 



section 111. We will see below, however, that the relic abundance will not increase when 



intermediate scales are considered. 



3.1 Inflation scenario 



Let us consider for example the SUSY hybrid inflation scenario studied in ref. |^ . There, 
the inflaton decay width can be computed with the result 

2 



V = — ('111 



(9) 



^For another alternative cosmological scenario with the potential of increasing the relic density see 
ref. [|l4|, where the decay of cosmic strings producing neutralinos is considered. 
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where (0) is the vacuum expectation value of the inflaton field, which is of the order of 
the unification scale, is the infiaton mass, and mj is the mass of the particle / that 
the infiaton decay to (in this case a right handed neutrino or sneutrino). Obviously, rrif 
should be smaller than the infiaton mass to allow for the decay ff. 

Now, using eqs. (||) and (|^) one obtains the following reheating temperature: 

r..-.r.u,a.v(l^) (^) {^f (^f. (.) 



In ref. WE\ it has been shown that an intermediate unification scale of the order of Mj ~ 
10^^ GeV is favoured by infiation. Then, recalling that the infiaton mass is constrained 
by ^(j> ^ M] /Mpianck, "we obtain ~ 10^ GeV. From Eg . (p!0D we find that Tr/^ ~ 1 
GeV, since now (0) ~ 10^^ GeV. This reheating temperature is lower than the typical 
freeze-out temperature Tp ~ m^/20. Notice that in the standard GUT scenario discussed 
in ref. [|T5[, one has ~ 10" GeV, (0) ~ 10^^ GeV, and therefore Trh ~ 10^ GeV. 

As mentioned above, a detailed analysis of the relic density with a low reheating 
temperature has been carried out in ref. [|Il|] by Giudice, Kolb and Riotto. They study 
two possible non-relativistic cases, depending on whether or not the dark-matter particles 
are in chemical equilibrium. In the first one they are never in equilibrium, either before or 
after reheating. In the second one the dark-matter particles reach chemical equilibrium, 
but then freeze out before the completion of the reheat process. These scenarios not only 
lead to different qualitative and quantitative predictions for the relic density, but also 
these predictions are quite different from the standard ones summarised in eq. (^. 

In the case of non-equilibrium production, the number density of neutralinos n^o is 
much smaller than n^o , thus the relevant Boltzmann equations can be approximated and 



solved. One gets |11 



^.----o^(|)=(^)-(J^)'ll2|;^(..^), ,n) 

where g is the number of degrees of freedom of the neutralino and T^, is the temperature 
at which most of the neutralino production takes place, it is given by ~ 4m^o/15. 
As we can see ^x^h'^ is proportional to the annihilation cross section, instead of being 
inversely proportional as in eq. (|^). This raises the hope that the relic abundance could be 
increased in scenarios with intermediate scales where generically it is low. Unfortunately, 
the assumption that n^o << n^o leads to a severe constraint on the annihilation cross 

section [|11|. Namely < and ap < Op, where and Op are of the order of 10~^^ 
GeV~^ and 10~^^ GeV"^, respectively, for Trh ~ 1 GeV. Since we are interested in large 
cross sections, of the order of 10~^ GeV~^, eq. (|1^) cannot be applied. 

With a large annihilation cross section (a^ > a.s or ap > dp), the neutralino reaches 



equilibrium before reheating as discussed in ref. [|rT|, and its relic density is given by 



2 .o...n-ii9y"(TRH) n^GeV-^ 



Q^oh' = 2.3 X iQ-ii ^- ^HH^ ^ (^2) 



Now the relic density is again inversely proportional to the annihilation cross section as in 
eq. (^). Moreover, it has a further suppression because of the low reheating temperature 
TjiH ~ 1 GeV, and as a consequence we expect a result even worse than the one obtained 
in the standard computation discussed below eq. (|^). Indeed, for a neutralino-nucleus 
cross section of the order of 10^^ GeV~^ we obtain Vt^oh? ^ 10~^. 



3.2 Modulus— decay scenario 

It has been assumed in the above computation that the relic abundance does not receive 



any contribution from other sources. However, as shown by Moroi and Randall [0, the 
production of neutralinos through moduli decay can modify those results. 

Let us recall first that moduli fields are present for example in string theory Since 
moduli acquire masses through SUSY breaking effects, these masses, rriip, are expected to 
be of the order of the gravitino mass, i.e. (9(100 — 1000 GeV). On the other hand, their 
couplings with the MSSM matter are suppressed by a high energy scale. Thus one can 
parameterise the moduli decay width as 



ml 



2nMf ' 



(13) 



where we denote with Mj the effective suppression scale. Since we are interested in the 
analysis of scenarios with intermediate unification scales, we will consider the case of 
Mj ^ 10^^"^"^ GeV. An explicit example where this situation arises is the case of type I 
string constructions. There, twisted moduli are present with interactions suppressed by 
the string scale. As discussed in ref. |TB| intermediate values for this scale can be obtained. 



and they have very interesting phenomenological implications ||18| , |T9|, ^, |T6[ . 



Using eqs. (^ and ([T3| ) one obtains the following reheating temperature: 

T,„ = 1 MeV X f^^li£!l^) f^n^ f (^\" . (14) 

where note that = 10.75 for T ~ 0{1 — 10) MeV. This reheating temperature is 
shown as a function of the modulus mass in Fig. |l| for different values of M/. The 
request that the modulus mass is larger than ~ 100 — 500 GeV in order to allow for 
kinematical decays into neutralinos of suitable mass m-j^o ~ 50—200 GeV, limits in practice 
the reheating temperature to be above ~ 3 GeV for the lowest scale on consideration 
M/ = 10^^ GeV. This has important consequences for the relic density computation. As 
discussed in the introduction of this section, we need a temperature smaller than the 
typical Tp ~ m^o/20 ~ 3 — 10 GeV, in order to increase the relic neutralino density. This 
implies that the scale Mj ~ 10^^ GeV is in the border of validity. On the other hand, for 
larger values we can obtain very easily interesting reheating temperatures. For example, 
for Mj = 10^^ GeV we have T^h ^ 0.3 GeV. For the highest scale with an interesting 
phenomenological value of the neutralino-nucleus cross section, in the case of universality 

•^See other examples of moduli fields, for instance in GUTs, in ref. |17|]. 
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Figure 1: The reheating temperature T^h as a function of the modulus mass m^. The 
six curves correspond, from top to bottom, to Mj = 10", 10^^ 10^^ 10^^, 10^^ 10^*^ GeV, 
respectively. The region bounded by the horizontal lines corresponds to a reheating 
temperature larger than 1 MeV, a value close to the the nucleosynthesis limit, and smaller 
than 1 GeV, a value close to the freeze-out limit, as explained in the text. 



and moderate tan/3 0, Mj = 10^^ GeV, the lowest value of the reheating temperature 
corresponds to Tr// ~ 6 MeV. Larger values of the scale, Mj > 10^^ GeV, producing also 
a large cross section, are possible in D-brane scenarios since non-universality in soft terms 
is generically present ||^. In this case the constraint Tjih ^ 1 MeV from nucleosynthesis 
can be translated into a constraint on the modulus mass > 140 GeV. 

When considering the decay of the modulus field producing neutralinos, the evolu- 
tion of the cosmological abundance of the latter becomes more complicated than in the 
usual thermal-production case reviewed in Section 2. Now one has to solve the coupled 
Boltzmann equations for the neutralino, the moduli field and the radiation ||13|, pO|: 



dn^o 

— ^ + 3Hn.o 
dt ^1 

dn^ 

+ 6Hns 



dt 

dprad 

dt 



+ 4.Hprad 



(15) 
(16) 

(17) 



where A^^o is the averaged number of neutrahnos produced in the decay of one modulus 
field. 



Let us discuss qualitatively the solution following the arguments used in ref. [13]. For 



a Trh higher than Tp the relic density will roughly reproduce the usual result given by 
eq. (^. However, for the interesting case for us when Trh is lower than Tp, neutrali- 
nos produced from modulus decay are never in chemical equilibrium, unlike the thermal 
production case reviewed in Section 2. As a consequence, its number density always de- 
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creases through pair annihilation. When the annihilation rate ((j|o"f )n^o drops below the 
expansion rate of the universe, if, the neutralino freezes out. Then the relic density can 



be estimated as \13 



il,.h- = , ^/'"g'-^ , (18) 



This result is valid when there is a large number of neutralinos produced by the modulus 
decay. When the number is insufficient, they do not annihilate and therefore all the 
neutralinos survive. The result in this case is given by 

^1 (27r745) pjs^ ^ ' 

The actual relic density is estimated ||13[ as the minimum of (0) and ([TOl) . 

Now we can apply the above equations to our case with intermediate scales. Using 
eqs. (|T3|) and (|14[), we can write expressions ([T8|) and ( pJf ) as 



^1 U. 5 X 1020 Gel// I ((T|rt^) M m,^ y V ^* / UOOGel^/ ' 



2_ ~ ,'1.2 X lO^o GeV\ / m,^ /l0.75\ ' / m^o 



From these equations we can see that even with a large annihilation cross section, (cr^o"f ) ~ 
10~^ GeV~2^ we are able to obtain the cosmologically interesting value ^^^o/i^ ~ 1. For 



example for Mj ~ 10 GeV we obtain it when A'^o ~ 1, using eq. (|20|) . In Fig. y we 
show in more detail these results solving numerically the Boltzmann eqs. (p!5|)-(|T^), for the 
large annihilation cross section introduced in eq. (|^). There, the contours of constant relic 
neutralino density ^^oh"^ as a function of and N^o are shown, for fixed values of M/ 
and iTL^o. In particular, we consider the cases Mj = 10^"^, 10^^, 10^^ GeV, with m^o = 100 
GeV. The corresponding reheating temperatures can be obtained from Fig. 1. Note that 
whereas many values of N^o correspond to a satisfactory relic density for Mj = 10^^ _ xg^^ 
GeV, for the case Mj = 10^^ GeV only a small range works. 

Let us finally remark that the numerical value of N^o is in general model dependent. 



This was discussed in the context of supergravity in ref. [T^l. In this particular case both 
values N^o ~ 1 and iV^o ~ 10^^ — 10~^ are plausible, depending on the characteristics of 
the supergravity theory under consideration. 
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Figure 2: relic neutralino density ^j^oh"^ contours as a function of and iV^o, for m^o = 
100 GeV and several possible values of the intermediate scale Mi — 10^^, 10^^, 10^'' GeV. 
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4 Conclusions 



Current dark matter experiments are sensitive to a large neutralino-nucleus cross section, 
cr^o_;v ^ 10~^-10~^ GeV~^. There are regions in the parameter space of SUSY scenarios 
with an intermediate unification scale, where these cross sections can be obtained. How- 
ever, in these regions, the standard computation of the relic abundance of neutralinos 
through thermal production in the early universe, would imply a small relic abundance, 

n^oh^ < 0.01. 

We have analysed some alternatives to solve this potential problem. Let us recall 
that in the standard computation one is tacitly assuming that the reheating temperature 
is much larger than the freeze-out temperature, Tp ~ 10 GeV, and originated in an 
inflationary process. We have shown however that, when intermediate scales are present, 
reheating temperatures as small as (9(1 — 1000 MeV) are possible. Unfortunately, although 
the result for the relic abundance is modifled, still Q^oh"^ is too small. On the other 
hand, when the above reheating temperatures are associated to the decay of moduli fields 
producing out of equilibrium neutralinos, the relic abundance increases considerably with 
respect to the standard production. Thus the neutralino becomes a good dark matter 
candidate with 0.1 < Qh"^ < 0.3, for intermediate scales M/ ~ 10^^ — lO^'' GeV, and 
moduh mass ~ 100 — 1000 GeV. 
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